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Abstract. We present a review of recent work on the mathematical 
aspects of the BCS gap equation, covering our results of [3] as well 
our recent joint work with Hamza and Solovej [8] and with Frank and 
Naboko [6], respectively. In addition, we mention some related new 
results. 



1. Introduction 

In this paper we shall describe our recent mathematical study [HI [6l [9] of 
one of the current hot topics in condensed matter physics, namely ultra cold 
fermionic gases consisting of neutral spin-^ atoms. The kinetic energy of 
these atoms is described by the non-relativistic Schrodinger operator, and 
their interaction by a pair potential XV with A being a coupling parameter. 
As experimentalists are nowadays able to vary the inter-atomic potentials, 
the form of in actual physical systems can be quite general; see the recent 
reviews in [5] and [1] . Our primary goal concerns the study of the superfluid 
phases of such systems. According to Bardeen, Cooper and Schrieffer [2] 
(BCS) the superfluid state is characterized by the existence of a non-trivial 
solution of the gap equation 



at some temperature T > 0, with E{p) = \J (p^ — /x)^ -|- | A(p)p. Here, 
/i > is the chemical potential and V{j)) = (27r)~^/^ Jj^s V{x)e~^P^dx de- 
notes the Fourier transform of V. The function A(p) is the order parameter 
and represents the wavefunction of the Cooper pairs. Despite the fact that 
the BCS equation ([1]) is highly non- linear, we shall show in Theorem [1] (see 
also [HI Thm 1]) that the existence of a non-trivial solution to ([1]) at some 
temperature T is equivalent to the fact that a certain linear operator, given 
in ([6]) below, has a negative eigenvalue. For T = this operator is given by 
\ — A — XV. This rather astonishing possibility of reducing a non-linear 
to a linear problem allows for a more thorough mathematical study. Using 
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spectral-theoretic methods, we are able to give a precise characterization of 
the class of potentials leading to a non-trivial solution for ([T]). In partic- 
ular, in Theorem [2] (see also [6l Thm 1]) we prove that for all interaction 
potentials that create a negative eigenvalue of the effective potential on the 
Fermi sphere (see ([9]) below; a sufficient condition for this property is that 
/ig3 V{x)dx < 0), there exists a critical temperature Tc{XV) > such that ([T]) 
has a non-trivial (i.e., not identically vanishing) solution for all T < Tc{XV), 
whereas there is no such solution for T > Tc{XV). Additionally, we shall 
determine in Theorem [2] the precise asymptotic behavior of Tc{XV) in the 
small coupling limit. We extend this result in Theorem [3] (see also [9, Thm 
1]) and give a derivation of the critical temperature Tc valid to second order 
Born approximation. More precisely, we shall show that 

T, = t,- e-/(2v^^M) (2) 

TT 

where 7 ^ 0.577 denotes Euler's constant, and where 6^ < is an effective 
scattering length. To first order in the Born approximation, 6^ is related 
to the scattering amplitude of particles with momenta on the Fermi sphere, 
but to second order the expression is more complicated. The precise formula 
is given in Eq. (jl4p below. For interaction potentials that decay fast enough 
at large distances, we shall show that 6^ reduces to the usual scattering 
length ao of the interaction potential in the low density limit, i.e., for small 
/i. Our formula thus represents a generalization of a well-known formula in 
the physics literature O [13] . 

In the case of zero temperature, the function E{p) in ([1]) describes an 
effective energy-momentum relation for quasi particles, and 

H := inf E{p) = inf ^7(^2 _ ^)2 + |A(p)|2 
p p 

is called the energy gap of the system. It is of major importance for applica- 
tions, such as the classification of different types of superfluids. In fact, H is 
the spectral gap of the corresponding second quantized BCS Hamiltonian. 
(See p] and [l2] or the appendix in [8j.) 

An important problem is the classification of potentials V for which S > 0. 
This questions turns out to be intimately related to the continuity of the 
momentum distribution 7(p), which will be introduced in the next section. 
In the normal (i.e., not superfluid) state, A = and 7 is a step function at 
T = 0, namely j{p) = 0{\p\ — ^/Jl). According to the picture presented in 
standard textbooks the appearance of a superfiuid phase softens this step 
function and 7(p) becomes continuous. We are going to prove in this paper 
that if y(x)|x| G L^/^ and Jy < then indeed both strict positivity of 
H > and continuity 7 hold. It remains an open problem to find examples 
of potentials such that the gap vanishes in cases where a superfluid phase 
occurs. 
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One of the difficulties involved in evaluating S is the potential non- 
uniqueness of the solution of the BCS gap equation. For interaction po- 
tentials that have nonpositive Fourier transform, however, we shall show 
that the BCS pair wavefunction is unique, and has zero angular momen- 
tum. In this case, we shall prove in Theorem [5] (see also [9, Thm. 2]) similar 
results for H as for the critical temperature. It turns out that, at least up 
to second order Born approximation, 

E = Tc— (3) 

in this case. This equality is valid for any density, i.e., for any value of 
the chemical potential /i. In particular, H has exactly the same exponential 
dependence on the interaction potential, described by 6^, as the critical 
temperature Tc. 

2. Preliminaries and main results 

We consider a gas of spin 1/2 fermions at temperature T > and chemical 
potential fJ. > 0, interacting via a local two-body interaction potential of the 
form 2XV{x). Here, A > is a coupling parameter, and the factor 2 is 
introduced for convenience. We assume that V is real-valued and has some 
mild regularity properties, namely V £ L^{R^) n L^/^(M^). In the BCS 
approximation, the system is described by the BCS functional J-t, derived 
by Leggett in his seminal paper [11], based on the original work of BCS [2]. 
The BCS functional J-t is related to the pressure of the system and is given 
by 



TT{l,a) = J {p'-^)j{p)dp + J \aix)\'Vix)dx-TS{j,a), (4) 
where the entropy S is 

5(7, «) = - / Tre. [Tip] log r(p)] dp, r{p) = (iM ^ 

The functions 7(p) and a{p) are interpreted as the momentum distribution 
and the Cooper pair wave function, respectively. The satisfy the matrix 
constraint < T{p) < 1 for all p G M^. In terms of the BCS functional 
the occurrence of superfluidity is described by minimizers with a ^ 0. We 
remark that in the case of the Hubbard-model this functional was studied 
inH. 

For an arbitrary temperature < T < oo the BCS gap equation, which 
is the Euler-Lagrange equation associated with the functional J^t, reads 

where E{p) = ^ (p^ — [i)^ + |A(p)p. The order parameter A is related to 
the expectation value of the Cooper pairs a via 2a(p) = A(p)/i?(p). We 
present in the following a thorough mathematical study of this equation. In 
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order to do so, we shall not attack the equation ^ directly, but exploit the 
fact that a is a critical point of the semi-bounded functional J-t- 

The key to our studies is the observation in [8j that the existence of 
a non-trivial solution to the non-linear equation ([5]) can be reduced to a 
linear criterion, which can be formulated as follows. 

THEOREM 1 ([g Theorem 1]). Let V G L^^, fi G R, and oo > T >0. 

Define 



,2 .(p'-^^yT + i 



Then the non-linear BCS equation ([5]) has a non-trivial solution if and only 
if the linear operator 

KT,f, + XV, (6) 
acting on L^(M^), has at least one negative eigenvalue. 



Hence we are able to relate a non-linear problem to a linear problem which 
is much easier to handle. The operator Kt,^ is understood as a multiplica- 
tion operator in momentum space. In the limit T ^ this operator reduces 
to \ - A- + XV. 

2.1. The critical temperature. Theorem [1] enables a precise definition of 
the critical temperature, by 

T,{XV) := inf{r |i^T,M + > 0}. (7) 

The symbol Kt,^{p) is point-wise monotone in T. This implies that for any 
potential V, there is a critical temperature < Tc{XV) < oo that separates 
two phases, a superfluid phase for < T < Tc{XV) from a normal phase for 
Tc{XV) <T<oo. Note that Tc{XV) = means that there is no superfluid 
phase for XV. Using the linear criterion ([7]) we can classify the potentials 
for which Tc{XV) > 0, and simultaneously we can evaluate the asymptotic 
behavior of Tc{XV) in the limit of small A. This can be done by spectral 
theoretical methods. Applying the Birman-Schwinger principle one observes 
that the critical temperature Tc can be characterized by the fact that the 
compact operator 

X{sgnV)\V\'/'K-}jV\'/' (8) 

has —1 as its lowest eigenvalue. This operator is singular for Tc 0, and 
the key observation is that its singular part is represented by the operator 
Aln(l/Tc)V;„ where : L\%) ^ L'^{%) is given by 

Here, denotes the 2-sphere with radius ^/Jl, and dto denotes Lebesgue 
measure on fi^. We note that the operator has appeared already earlier 
in the literature [H [10] . 

Our analysis here is somewhat similar in spirit to the one concerning 
the lowest eigenvalue of the Schrodinger operator p'^ + XV in two space 



SPECTRAL PROPERTIES OF THE BCS GAP EQUATION OF SUPERFLUIDITY 5 



dimensions [H]. This latter case is considerably simpler, however, as has 
a unique minimum at p = 0, whereas Kt,^{p) takes its minimal value on 
the Fermi sphere = /x, meaning that its minimum is highly degenerate. 
Hence, in our case, the problem is reduced to analyzing a map from the 
functions on the Fermi sphere (of radius to itself. Let us denote the 
lowest eigenvalue of as 

eM(^) := infspecV^. 

Whenever this eigenvalue is negative then the critical temperature is non 
zero for all A > 0, and we can evaluate its asymptotics. Moreover, the 
converse is "almost" true: 

THEOREM 2 ([6, Theorem 1]). Let V G L3/2(m3) p^i(]^3) real-valued, 
and let A > 0. 

(i) Assume that e^{V) < 0. Then Tc{XV) is non-zero for all A > 0, and 

lim A In t^ttttt^ = tttt • (10) 

A^o T,{XV) e^{V) ^ ^ 

(ii) Assume thate^{V) = 0. IfTc{XV) is non-zero, then\n{fi/Tc{XV)) > 
cA~^ for some c > and small X. 

(iii) // there exists an e > such that e^(y — e\V\) = 0, then Tc{XV) = 
for small enough X. 

As we see, the occurrence of superfluidity as well as the asymptotic be- 
havior of Tc{XV) is governed by e^{V). A sufficient condition for e^{V) to 
be negative is J ^ < 0. But one can easily find other examples. Eq. (jlOp 
shows that the critical temperature behaves like Tc{XV) ~ fie^^^'^'^'^^^^h In 
other words it is exponentially small in the coupling. 

In the following, we shall derive the second order correction, i.e., we will 
compute the constant in front of the exponentially small term in Tc- For 
this purpose, we define an operator on L^(f]^) via its quadratic form 

{u\w,\u) = d\p\ (^^y^ [^^ {Wip)? - \V{^P/\P\)?) 

l,dn\^(^p/\p\)f^ . (11) 

Here, ip{p) = (27r)-3/2 J^^ y(p _ q)u{q)dw{q), and {\p\,n) e M+ x de- 
note spherical coordinates for p £ M^. We note that since V £ L^(]R^), 
Jg2 dQ |(^(p)p is Lipschitz continuous in \p\ for any u £ L^(M'^), and hence 
the radial integration is well-defined, even in the vicinity of p^ = fi. In fact 
the operator can be shown to be Hilbert-Schmidt class, see ^9^ Section 
3]. 

For A > 0, let 
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and let 6^ (A) denote its ground state energy, 

6^(A) = inf specS^ . (13) 

We note that if < 0, then also 6^ (A) < for small A. In fact, if the 
eigenfunction corresponding to the lowest eigenvalue of is unique and 
equals u G Lp'iTt^), then 

h,{X) = {u\BM + O(A^) = - A^^^i^i^ + O(A^) . (14) 

In the degenerate case, this formula holds if one chooses u to be the eigen- 
function of that yields the largest value (M|W^|ii) among all such (nor- 
malized) eigenf unctions. 

With the aid of 6^ (A), we can now recover the next order of the critical 
temperature for small A. 

THEOREM 3 ([1 Theorem 1]). Let V G L'^{W^)rM^I'^{m?) and let fi>0. 
Assume that e^ = inf spec V^^ < 0, and let b^{X) be defined in il3\). Then the 
critical temperature Tc for the BCS equation is strictly positive and satisfies 

lim fin + ^^TTT^ = 2 - 7 - ln(8/7r) . (15) 

Here, 7 ~ 0.577 denotes Euler's constant. 
The Theorem says that, for small A, 

TT 

Note that 6^ (A) can be interpreted as a (renormalized) effective scattering 
length of 2XV{x) (in second order Born approximation) for particles with 
momenta on the Fermi sphere. In fact, if V is radial and J^s V{x)dx < 0, 
it is not difficult to see that for small enough fi the (unique) eigenfunction 
corresponding to the lowest eigenvalue of is the constant function 
u{p) = (47r^)-i/2. (See [6, Section 2.1].) For this u, we have 

lim{u\BJu) = (A/47r) [ V{x)dx - {X/Airf [ Y^^^^Mfixdy = ao(A) . 
M^O J^a J^6 \x - y\ 

Here, ao(A) equals the scattering length of 2Ay in second order Born ap- 
proximation. Assuming additionally that y(x)|x| G and bearing in mind 
that 6^(A) = {u\B^\u) +0{\^) for small enough fi, we can, in fact, estimate 
the difference between 6^ (A) and ao(A). Namely we prove in ^ Proposition 
1] that 

hm —— — — = . 

i^-*o ^ \{u\B^\u) ao[X)J 

This yields the approximation 

~ ^^£Zl!g7r/(2yAlao(A)) 

vr 
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in the limit of small A and small /i. This expression is well-known in the 
physics literature [3 [13]. We point out, however, that our formula (fT6]l is 
much more general since it holds for any value of /i > 0. 

2.2. Energy Gap at Zero Temperature. Consider now the zero tem- 
perature case T = 0. In this case, it is natural to formulate a functional 
depending only on a instead of 7 and a. In fact, for T = the optimal 
choice of jip) in J^t for given d{p) is clearly 

_ M(l + Vl-4|a(p)P ) for p2 < ^ 
1 i(l-Vl-4|a(p)P) forp2>^ • 

Subtracting an unimportant constant, this leads to the zero temperature 
BCS functional 

•^o(«) = ^/" (1- Vl-4|a(p)|2)dp+A / y(x)|a(x)|2dx. (18) 

The variational equation satisfied by a minimizer of (jlSp is then 

with A(p) = 2£'(p)a(p). This is simply the BCS equation ([2D at T = 0. For 
a solution A, the energy gap H is defined as 



E = inf ^(p) = inf v/(p2 -//)2 + |A(p)|2 . (20) 
p p 

It has the interpretation of an energy gap in the corresponding second- 
quantized BCS Hamiltonian (see, e.g., [12] or the appendix in [8j.) 

A priori, the fact that the order parameter A is non vanishing does not 
imply that H > 0. Strict positivity of H turns out to be related to the 
continuity of the corresponding 7 in (|17p . In fact, we are going to prove in 
Lemma [1] that if V decays fast enoug h, i.e., V{x)\x\ E L^/^{R^), the two 
properties, H > and 7(p) continuous, are equivalent. Both properties hold 
true under the assumption that J V < 0: 

THEOREM 4. Let V G L^/^ n L^, with V{x)\x\ £ L^/^{R^) and JV = 
(27r)'^/2y(Q) < j^^i a be a minimizer of the BCS functional. Then H 
defined in i20\) is strictly positive, and the corresponding momentum distri- 
bution 7 in (11) is continuous. 

One of the difficulties involved in evaluating H is the potential non- 
uniqueness of minimizers of (jlSp . and hence non-uniqueness of solutions 
of the BCS gap equation (I19p . The gap H may depend on the choice of A 
in this case. For potentials V with non-positive Fourier transform, however, 
we can prove the uniqueness of A and, in addition, we are able to derive the 
precise asymptotic of S as A ^ 0. 

In the following we will restrict our attention to radial potentials V with 
non-positive Fourier transform. We also assume that V'(O) = (27r)~^/2 j y(^x)dx < 
0. It is easy to see that = inf spec < in this case, and that the 
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(unique) eigenfunction corresponding to this lowest eigenvalue of is the 
constant function. 

In particular we have the following asymptotic behavior of the energy gap 
S as A — > 0. 

THEOREM 5 ([9, Theorem 2]). Assume that V G L^{M.^) D L^^^iR^) is 
radial, with V{p) < and V{0) < 0. Then there is a unique minimizer (up 
to a constant phase) of the BCS functional M^) atT = 0. The corresponding 
energy gap, H = infp \/ (p^ — /i)^ + |A(p)p , is strictly positive, and satisfies 

lim ( In f + ^^vTT ) = 2 - ln(8) . (21) 

Here, b^{X) be defined in U3\) . 

The Theorem says that, for small A, 

In particular, in combination with Theorem [3l we obtain the universal ratio 

lim = — « 1.7639. 

A-.0 Tc e-y 

That is, the ratio of the energy gap H and the critical temperature Tc tends to 
a universal constant as A — > 0, independently of V and fi. This property has 
been observed before for the original BCS model with rank one interaction 
[21 [T2] , and in the low density limit for more general interactions [7] under 
additional assumptions. Our analysis shows that it is valid in full generality 
at small coupling A ^ 1. 

3. Sketch of the proof of Theorem [T] 

The backbone of our analysis is the linear criterion in Theorem [TJ As a 
first step towards its proof, one has to prove that the functional J-'xi'y, a) in 
([H) attains a minimum on the set 

V = {(7,a)|7 G L^(M^(l+p2)dp),a g H^{R^),0 < 7 < 1, |a|2 < 7(1-7)}. 

This can be done by proving lower semi-continuity of J-t on T>. See [U 
Prop. 1] for details. Theorem [T] is then a direct consequence of the equiva- 
lence of the following three statements [8, Theorem 1]: 

(i) The normal state (70, 0), with 70 = [e'^^^^'^^^-^-l-l]"^ being the Fermi- 
Dirac distribution, is unstable under pair formation, i.e., 

inf Trh, a) < J^rilo, 0) . 

(ii) There exists a pair (7, a) £ V, with a 7^ 0, such that 

A(p) = -^^^Mp) (22) 
satisfies the BCS gap equation ([5]). 
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(iii) The linear operator Kt,^ + V has at least one negative eigenvalue. 

The proof of the equivalence of these three statement consists of the fol- 
lowing steps. First, it is straightforward to show that (i) =^ (ii). By eval- 
uating the stationary equations in both variables, 7 and a, one shows that 
the combination (j22p satisfies the BCS equation 

To show that (iii) (i), first note that (70,0) is the minimizer of Tt in 
the case V = Q. Consequently ^J^T{lOi'tg)\t=o = for general g. Moreover, 
a simple calculation shows that 

-^Hlo,tg)t=o = 2{g\KT,^ + XV\g) . 

If KT,fj.+\V has a negative eigenvalue, we thus see that J^{'~fo,tg) < .Ft(7o, 0) 
for small t and an appropriate choice of g. 

The hardest part in showing the equivalence of the three statements is to 
show that (ii) (iii). Given a pair (7,0) such that the corresponding A 
in (f22|) satisfies the BCS equation ([5]), we note that if a = m{p)dt{p) and 
7(p) = 1/2 + m{p){'y{p) — 1/2), the pair (7, a) yields the same A and hence 
also satisfies ([5]). Moreover, with the choice 

tanh^ 

m{p) = - 



2 



7(p) 2E{p) 



(where E{p) = \J (jp' — fi)"^ + | A(p)P), the new pair (7,0;) satisfies addition- 
ally 

2Eip) p^-fi 



tanh^^ ^-7(P) 



(23) 



^jvip- q)a{q)dq = - ^^Hp) ■ (24) 



2 

Note that in the case V = 0, i.e., A = 0, the equation ([23|) reduces to 

2-70 

Using this fact, together with (I24p . we thus obtain 



{a\KT,^, + \V\a) 



1 9 
p — fi p — ^ 



1-70 \-l 



a ) . (25) 



Using the definition of E{p) and the strict monotonicity of the function 
X I— > x/tanh ^ for x > 0, we infer from (j23p that 



- ^^ < P^ ~ ^^ 
I - 70 ~ \ 



with strict the inequality on the set where A 7^ 0. Consequently, the expres- 
sion (p5|) is strictly negative. Hence i^^r.^i + has a negative eigenvalue. 
This shows that (ii) implies (iii). 



10 



CHRISTIAN HAINZL AND ROBERT SEIRINGER 



4. Proof of Theorems [2] and [3] 

For a (not necessarily sign-definite) potential V{x) let us use the notation 

y(x)i/2 ^ {sgnV{x))\V{x)\^/^ . 

Prom our definition of the critical temperature Tc it follows immediately 
that for T = Tc the operator Kt,^ + XV has and eigenvalue and no neg- 
ative eigenvalue. If tp is the corresponding eigenvector, one can rewrite the 
eigenvalue equation in the form 

Multiplying this equation by V^^'^{x), one obtains an eigenvalue equation 
for (f = This argument works in both directions and is called the 

Birman-Schwinger principle (see [U Lemma 1]). In particular it tells us that 
the critical temperature Tc is determined by the fact that for this value of 
T the smallest eigenvalue of 

Bt = XV^/^K:^^^\V\^/^ (26) 

equals —1. Note that although Bt is not self-adjoint, it has real spectrum. 

Let ^ : L^(M^) Lp'iQ.^) denote the (bounded) operator which maps 
■0 G L^(M^) to the Pourier transform of restricted to the sphere fi^. Since 
V eL^{M?), multiplication by \ V\^/'^ is a bounded operator from L^(M'^) to 
L1(M3), and hence ^\V\^/'^ is a bounded operator from L^(]R^) to L'^{Q^). 
Let 

r-.m = {i^ I (if^^ - ^) * . »} . 

and let 

MT = K^l-m^{T)r^- (27) 
As in [61 Lemma 2] one can show that y^/2^^|y|i/2 

is a Hilbert-Schmidt 

operator on 

L2(m3)^ and its Hilbert Schmidt norm is bounded uniformly in 

T. In particular, the singular part of Bt as T — > is entirely determined 
by yV2^*;j|y|i/2. 

Since V^^'^Mt\V\^^'^ is uniformly bounded, we can choose A small enough 
such that 1 + Ayi/2Mr|y|^/^ is invertible, and we can then write 1 + Bt as 

1 + Bt = 1 + Ayi/2 (m^(T)d*d + Mt) \V\^/^ (28) 
Then Bt having an eigenvalue —1 is equivalent to 



(29) 
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having an eigenvalue —1. The operator in (f29l) is isospectral to the selfadjoint 
operator 

acting on L'^{Q^). 

At T = Tc, —1 is the smallest eigenvalue of Bt, hence (p9|) and ([30|) have 
an eigenvalue —1 for this value of T. Moreover, we can conclude that —1 is 
actually the smallest eigenvalue of (j29p and (j30p in this case. For, if there 
were an eigenvalue less then —1, we could increase T and, by continuity, find 
some T > Tc for which there is an eigenvalue —1. Using (j28p . this would 
contradict the fact that Bx has no eigenvalue —1 for T > T^. 

Consequently, the equation for the critical temperature can be written as 

Arn,(T.)infspecg|y|V^ ^^^^^^,^^^^^^^^^, yV2y = _i. (31) 

This equation is the starting point for the proof of Theorems [2] and [3l 
Proof of Theorem Up to first order in A the equation (j3ip reads 

Xnif^iTc) mispecdiV - XV Mt^V + 0{X'^)]r = -1 , (32) 

where the error term O(A^) is uniformly bounded in T^. Note that ^V^* = 
■y/JIVfj, defined in ([9]). Assume now that = inf specV^ is strictly negative. 
Since V^^'^Mt^V^^'^ is uniformly bounded, it follows immediately that 

lim Xmn(Tc) = = — — . 

A^o ' mfspecS'l^^* ^/J^e^ 

Together with the asymptotic behavior nif^iT) ~ ln(/i/T) as T — > 0, 

this implies the leading order behavior of ln(/i/rc) as A — > and proves the 
statement in (i). 

In order to see {ii) it suffices to realize that, in the case ^V^* > 0, Eq. ([32]) 
yields rriT^ > const / A^ . 

The statement (iii) is a consequence of the fact that 

^'^'^^ \ + AyV2M,jy|V2 ^^^^^* > m - const x\v\]r > 0, 

for A small enough. We refer to [B] for details. □ 

Proof of Theorem 0. To obtain the next order, we use Eq. (j32p and employ 
first order perturbation theory. Since ^V^* is compact and inf spec ^V^* < 
by assumption, first order perturbation theory implies that 

""^^^"^ " X{u\Wr\u) - X^{u\dVMT^Vd*\u) + 0(A3) ' ^^^^ 

where u is the (normalized) eigenfunction corresponding to the lowest eigen- 
value of ^V^* . (In case of degeneracy, one has to the choose the u that 
minimizes the A^ term in the denominator of (j33p among all such eigenfunc- 
tions.) 
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Eq. ()33p is an implicit equation for Tc. Since ^VMtV^* is uniformly 
bounded and Tc ^ as A —> 0, we have to evaluate the limit of {u\^VMtV'S*\u) 
as T ^ 0. To this aim, let = V^*u. Then 



{u\dVMTVd*\u) 

1 



j^^^^^^mp)rdp-m,{T)J^ \0{prduj{p) (34) 

-^-i^ [mpt - m^p/\p\r] + ^i^(v^p/n)p) dp. 

Recall that Kt,^{p) converges to — //| as T — > 0. Using the Lipschitz 
continuity of the spherical average of \ip{p)\^ (see P Eq. (29)]) it is easy to 
see that 

liin {u\WMtVT\u) = {u\Wf,\u) , (35) 

with yy^ defined in ([IT]) . In particular, combining ([33|) and ([35]) . we have 
thus shown that 

lim ( m^{Tc) + — . -^—^ 1=0. (36) 

A-o\ mfspec (A^V^-A2W^) / 



The statement follows by using the asymptotic behavior ([9l, Lemma 1]) 

m,{T) = -L (^in I + ^ - 2 + In ^ + o(l)) (37) 
in the limit of small T, where 7 ~ 0.5772 is Euler's constant. □ 



5. Proof of Theorems [H and [5] 

5.1. Sufficient condition for H > 0. If e^{V) < we know that the BCS 
equation has a solution, meaning the system shows a superfluid phase 
for r = 0. This is not sufficient, however, to guarantee the existence of 
a positive gap H > nor the continuity of the momentum distribution 7. 
Unlike the case of the critical temperature, we lack a linear criterion which 
allows a precise characterization of potentials V giving rise to a strictly posi- 
tive gap. We are, however, able to derive sufficient conditions, namely a fast 
enough decay of V. Under such assumptions one can show the equivalence 
of the positivity of H and the continuity of 7. Both hold true if additionally 
J V < 0. It remains an open problem to find examples for V such that 
< but S = 0. 

Lemma 1. Assume that V G and that V[x)\x\ G L^I^{R^). Then 

E > if and only if 7 is continuous. 

Proof. It is easy to deduce |8j from the BCS equation ([19]) that a is in 
C*^(]R^). Because of ([TT]) the continuity of 7 is equivalent to the fact that 
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Idl = 1/4 on the Fermi 17 



^. From the relation A(p) 
1 1 
4 ^ ■ 



\a{p)\' 



(p2 



|A{p)| 



^+1 



2E{p)a{p) one obtains 
(38) 



and we can conclude that \a\'^ = 1/4 on the Fermi surface if and only if 
A{p) does not vanish on fi^. Namely, suppose that A vanishes at some p' 
on the Fermi surface. Since a £ H^{R^) we see that a G L2(M3) n ^^(M^) 
and hence, together with y(j;)|j;| E L^^^, Holder's inequality implies that 
A(x)|x| = V{x)a{x)\x\ e L^(M'^). We thus infer that A(p) is Lipschitz 
continuous, meaning that A{p) cannot decay slower to than linear. Hence 

there is a 5 such that lim^^p/ ^^^^p > ^ and |a(p')p < 1 , < t- □ 



Proof of Theorem Let a be a global minimizer of the BCS functional J^q 
Then for any g £ Cq°(M^) such that |a + e^| < 1/2 for e small enough, 

jj='{a + eg) > 0. 



de2- 

A straightforward calculation yields 



■J^{a + eg) 



2{g\E{-iV)+XV\g) + 8 



\p^-M&m 



(39) 



(40) 



de^ ' ' ' J [l-4|a|2]3/2 

Assume now that H = 0. This means that A has to vanish at some point 
p' G 0^. Then there has to be an open neighborhood on on which 
A vanishes. In fact, according to the argument in the proof of Lemma [J 
(Eq. (jSSp and Lipschitz continuity of A) there is a neighborhood Ms{p') C 
in the vicinity of p' where |ap < 1/4—5 for some 5 > 0, and hence A vanishes 
on Ns{p') n il^. Note that A cannot vanish at one point on the Fermi 
surface since otherwise \a\ = 1/2 except on one point, which contradicts the 
continuity of a. 

We shall now construct an appropriate trial sequence gn, essentially sup- 
ported in Ms, such that 



lim 

n— >oo 



and 



{gn\E{-iV)\gn)+S 



lim {gn\y\gn) 



|p^-H[3^(«gn)]^ 
[l-4|d|2]3/2 

V{x)dx < 0. 







(41) 



(42) 



This gives a contradiction to ([39|) . 

For the construction of let V'n S L'^{Q^) be supported in Msip') n 17^ 
such that i^ni^s) — > 5{s — p') as n ^ oo. Choose also /„ G L'^{^+,t^dt) 
such that fn{t) ^(\//" — t), and let gn{p) = ipn{s) fn{\p\) ■ Observe that 
on M{p'), E{p) = \p^ — /i|/-y/l — 4|q;(p)P < cjp^ — ^| for some constant c, 
and thus grows linearly in \p\ close to -y//!. Hence one easily sees that the 
problem here is equivalent to the existence of a negative eigenvalue of the 
relativistic operator \p\ +V in one dimension. Using the Birman-Schwinger 
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principle, it is easy to see that the latter always has a negative eigenvalue if 
JV<0. □ 

5.2. Proof of Theorem [5l The energy gap of the system at zero temper- 
ature, S = mipE{p), with 



Eip) = |p2 _ ;.|/Vl-4|a(p)|2 = V|p2_^|2 + |A(p)|2, 

depends on the behavior of |A(p)| on the Fermi sphere. The function A is 
not unique, in general and need not be radial even in case V is radial. 

Under the assumption that V is non-positive and V^(0) < 0, we shall argue 
in the following that the minimizer of the BCS functional (jlSp at T = is 
unique [U Lemma 3]. If, in addition, V is radial, this necessarily implies 
that also the minimizer has to be radial. Since < 0, 



a{p)V{p- q)a{q)dpdq > \a{pWiP - Q)\oiiQ)\ dpdq . (43) 



Hence, if a{p) is a minimizer of J^q, (fTSj) . so is |a(p)|. 

Assume now there are two different minimizers f ^ g, both with non- 
negative Fourier transform. Since t — > 1 — ^/l — At is strictly convex for 
< t < 1/2 we see that i/j = -^f + i^d^ satisfies 

^o(V') < \W) + \^o{9) ■ 

This is a contradiction to f^g being distinct minimizers, and hence f = g- 
In particular, the absolute value of a minimizer has to be unique. If a 
is the unique non-negative minimizer, then one easily sees from the BCS 
equation (using j V < Q) that a is, in fact, strictly positive. Hence any 
minimizer is non- vanishing. But (j43p is strict for non-vanishing functions, 
unless a{p) = |d(p)|e*'' for some constant k G M. 

To summarize, we have just argued that for F < 0, V^(0) < and V 
radial, the solution of the BCS equation is unique, up to a constant phase, 
and it is radially symmetric. This will enable us to apply the same methods 
as we used for the critical temperature in order to derive the asymptotic 
behavior of H. 

The variational equation (jl9|) for the minimizer of J-q can be rewritten in 
terms of a as 

{E{-iV) + \V{x))a{x) = Q. (44) 

That is, a is an eigenfunction of the pseudo differential operator E{—iV) + 
\V{x), with zero eigenvalue. Since ^ < and a(p) is non-negative we can 
even conclude that a has to be the ground state. 

Similarly to the proof of Theorem [31 we can now employ the Birman- 
Schwinger principle to conclude from (I44p that (px = V^^^a satisfies the 
eigenvalue equation 

Al-^/^^==^==|V^|^/^<^A = -</'A . (45) 
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Moreover, there are no eigenvalues smaller than —1 of the operator on the 
left side of g5]). 
Let 

mJA) = max J / ( ^ - dp , o \ . (46) 

Similarly to ([27|) . we split the operator in ([15|) as 

V'/'—^\V\'/' = m,iA)V'/'rd\V\'/' + V'/HIa\V\'/' . 
hi[—i\/) 

Again one shows that V^/^MaI^I^^^ is bounded in Hilbert-Schmidt norm, 
independently of A. Moreover, as in the proof of Theorem [3] (cf. Eqs. (j28p - 
pm i. the fact that the lowest eigenvalue of \V'^/'^E{-iV)~'^\V\^/^ is -1 is, 
for small enough A, equivalent to the fact that the selfadjoint operator on 

has —1 as its smallest eigenvalue. This implies that limx^o Xfh^( A) = 
— l/(y^e^) and hence, in particular, m^(A) ~ as A — > 0. The unique 
eigenfunction corresponding to the lowest eigenvalue < of is, in fact, 
a positive function, and because of radial symmetry of V it is actually the 
constant function u{p) = (47r//)~^/^. 

We now give a precise characterization of A{p) for small A. 

Lemma 2. Let V £ n l?!'^ he radial, with V < and V{0) < 0, and 
let A be given in U9\) . with a the unique minimizer of the BCS functional 
(d). Then 



A{p) = -/(A) ( [ 



Vip-q)dcj{q)+X7]xip)] (48) 



for some positive function /(A), with ||r?Al|L°°(R3) bounded independently of 
A. 

Proof. Because of (|45]) . ^?|T^|^/^</>a is the eigenfunction of ([47|) corresponding 
to the lowest eigenvalue —1. Note that because of radial symmetry, the 
constant function u{p) = (47r/i)~^/^ is an eigenfunction of (I47p . For small 
enough A it has to be eigenfunction corresponding to the lowest eigenvalue 
(since it is the unique ground state of the compact operator ^V^*). We 
conclude that 

for some normalization constant /(A). Note that ||^a||2 uniformly bounded 
for small A, since both V^/'^Ma\V\'^/'^ and F^s^* are bounded operators. 
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From ()44|) and the definition (px = V^^'^a we know tliat 

A(p) = 2E{p)a{p) = -2XVa{p) = -2X\Vy/^x{p) ■ 
In combination with ()49p this imphes that 

A{p) = -2A/(A) (vr^{p) + X^Txip)) , 

with r,x = \V\^/^^x. With llrfAlloo < (27r)-3/2||^^||^ < (27r)"3/2||y||i||^;,||2 by 
Schwarz's inequality, we arrive at the statement of the Lemma. □ 

With the aid of Lemma [2] and Lipschitz continuity of V{p — q) dw{q) 
(which follows from V € L^(M^)) it is not difficult to see that 



(50) 



as A — > 0. From Eq. ()47p we now conclude that 

"^^^^^ " \{u\Wr\u) - X^{u\WM/^V^*\u) + 0(A3) ' ^^'^^ 

where u{p) = (47r/i)~^/2 is the normalized constant function on the sphere 
O^. Moreover, with (p = V^*u, 

{u\WM^Vr\u) = j^^ dp - m^(A) ^ \^{^p/\p\)\'' du{p) 



L {e{p) 



^ [mp)\' - mvj^p/\p\)\'] + \mvj^p/\p\)\' ) dp. 



E{p) p' 
Using Lemma [2] and the fact that lim;^^o fW = Oi conclude that 

l[m{u\dVMAVd*\u) = {u\W^\u) , (52) 

A — *0 

with defined in pip . (Compare with Eqs. (j34p and (|35p .) In combination 
with (I50p and (I51|) and the definition of in (I12p , this proves that 

lim ( In ( -TT^ ) + — ^ , , 1 = 2 - ln(8) . 

The same holds true with {u\Bf^\u) replaced by bf^{X) = infspeC;B^, since 
under our assumptions on V the two quantities differ only by terms of order 
A3. 

Now, by the definition of the energy gap S in (120p . H < A(y^). Moreover, 

H > min |A(p)| , 

|p2-H<H' 

from which it easily follows that H > A(y7^)(l — o(l)), using Lemma[2l This 
proves Theorem [5l 
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